Excitability mediated by localized structures 
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We find and characterize an excitability regime mediated by localized structures in a dissipative 
nonlinear optical cavity. The scenario is that stable localized structures exhibit a Hopf bifurcation 
to self-pulsating behavior, that is followed by the destruction of the oscillation in a saddle-loop 
bifurcation. Beyond this point there is a regime of excitable localized structures under the applica- 
tion of suitable perturbations. Excitability emerges from the spatial dependence since the system 
does not exhibit any excitable behavior locally. We show that the whole scenario is organized by a 
Takens-Bogdanov codimension-2 bifurcation point. 

PACS numbers: 42.65.Sf, 05.45.-a, 89.75.Fb 



Localized structures (LS) in dissipative media have 
been found in a variety of systems, such as chemical re- 
actions, gas discharges or fluids 1]. They are also found 
in optical cavities, due to the interplay between diffrac- 
tion, nonlinearity, driving, and dissipation @, 0- These 
structures, also known as cavity solitons have to be dis- 
tinguished from conservative solitons found for example 
in propagation in fibers, for which there is a continuous 
family of solutions depending on their energy. Instead, 
dissipative LS are unique once the parameters of the sys- 
tem have been fixed. This fact makes this structures po- 
tentially useful in optical (i.e., fast and spatially dense) 
storage and processing of information 0, S • 

Localized structures may develop a number of insta- 
bilities like start moving, breathing or oscillating. In the 
latter case, they would oscillate in time while remaining 
stationary in space, like the oscillons found in a vibrated 
layer of sand 6]. The occurrence of these oscillons in 
autonomous systems has been reported both in optical 
0, E3, and chemical systems [Tof . 

In the present work we report on a route in which au- 
tonomous oscillating localized structures are destroyed, 
leading to a excitability regime. Excitability is a concept 
arising originally from biology (e.g. neuroscience) , and 
it has been found in a variety of systems ^]|, includ- 
ing optical systems 

nana. 

Typically a system is said 
to be excitable if while it sits at an stable fixed point, 
perturbations beyond a certain threshold induce a large 
response before coming back to the rest state. In ad- 
dition, excitability is also characterized by the existence 
of a refractory time during which no further excitation 
is possible. In phase space 0, excitability occurs 
for parameter regimes where a stable fixed point is close 
to a bifurcation in which an oscillation is created. Pos- 
sibly the best known example of an excitable system is 
FitzHugh-Nagumo model, close to the Hopf bifurcation, 
although one may also find excitable behavior mediated 
by a saddle point, namely either in the form of an An- 
dronov (or saddle-node on the invariant circle) bifurca- 
tion or a saddle-loop (or homoclinic) bifurcation. These 



three possible scenarios are the simplest possible, occur- 
ring in systems that, minimally, can be characterized by 
two phase variables. 

The concept of excitability has been extended to sys- 
tems with spatial dependence by coupling several or 
many zero-dimensional excitable systems Here we 

consider a different situation: a system that without spa- 
tial dependence does not show excitable behavior while 
the localized structures that appear in the spatially de- 
pendent systems do. 

A prototype model in nonlinear optics displaying the 
formation of LS is the one introduced by Lugiato and 
Lefever for an optical cavity filled with a Kerr medium 
[lq . In the mean field approximation the dynamics of 
the scaled slowly varying amplitude of the complex field 
E(x, t) is given by 



dE 



-(1 + iO)E + iV 2 E + E + i\E\ 2 E, 



(1) 



where V 2 = d 2 jdx 2 + d 2 /dy 2 . The first term describes 
cavity losses, E is a homogeneous (plane wave) driv- 
ing field, and the cavity detuning with respect to Eq. 
Eq. JU has been normalized by the cavity decay rate, 
and it has a homogeneous steady state solution E s given 
implicitly by E s = E /(l + i(0 - J s )), where I s = \E S \ 2 . 
We use the intra-cavity background intensity I s together 
with as convenient control parameters. The homoge- 
neous solution has a instability leading to the formation 
of a hexagonal pattern at I s = 1. The bifurcation starts 
subcritically and t he p attern coexists with the homoge- 
neous solution [Tfl, 0j| , so that a stable-unstable pair of 
LS is created in a saddle- node bifurcation 5]. The un- 
stable branch (the so-called lower or middle branch) has 
a single unstable direction. The region of existence of 
these LS, also known as Kerr cavity solitons, has been 
characterized in 0,07 an d is partially shown in Fig. E] 
Increasing the LS undergoes a supercritical Hopf bi- 
furcation and starts to oscillate autonomously 0, EH • 
For one spatial dimension, Eq. (QJ also exhibits LS in 
the appropriate parameter regime, but these structures 
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FIG. 1: a) Phase diagram: I s vs. LS are stable in the 
shaded region and oscillate in the dashed one. Lines indi- 
cate bifurcations: saddle- node (solid), Hopf (dot-dashed), and 
sadle-loop (dashed), b) Distance between the saddle- node and 
Hopf lines. 

do not undergo any Hopf bifurcation. 

As the control parameter is further increased part 
of the limit cycle moves closer and closer to the lower- 
branch LS, which is a saddle point in phase space, as 
illustrated in Fig. [21 On the left column we plot the time 
evolution of the LS maximum obtained from numerical 
integration of Eq. JQ , the dashed line shows for compar- 
ison the maximum of the lower-branch LS. On the right 
column we sketch the evolution on phase space projected 
on two variables. At a certain critical value a global bi- 
furcation takes place: the cycle touches the lower-branch 
LS and becomes a homoclinic orbit (Fig. (2t). This is 
an infinite-period bifurcation called saddle-loop or homo- 
clinic bifurcation [l9|. For > # c , the saddle connection 
breaks and the loop is destroyed (Fig. EJi). After fol- 
lowing a trajectory in phase space close to the previous 
loop the LS approaches the saddle point (dashed line) 
where the evolution is dominated by its slow stable man- 
ifold (see the long plateau between t = 15 and t = 60 in 
Fig. Eli). Finally the LS decays to the homogeneous so- 
lution (dotted line). For larger values of 6 the trajectory 
moves away from the saddle and, therefore, the decay to 
the homogeneous solutions takes place in shorter times. 

The saddle-loop bifurcation has a characteristic scaling 
law that governs the period T of the limit cycle as the 
bifurcation is approached. Close to the critical point the 
system spends most of the time close to the unstable 
LS (saddle). T can be then estimated by the linearized 
dynamics around the saddle 

Toc--^ln(0 c -0), (2) 

where Ai is the unstable eigenvalue of the saddle point. 
We are now going to show that this scaling law is verified 
in our system. Fig. 03 shows the period of the LS limit cy- 
cle as a function of the control parameter 0. As expected, 
the period of the limit cycle diverges logarithmically as 
the bifurcation is approached. We then evaluate Ai with 
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FIG. 2: Left: LS maximum intensity as a function of time 
for increasing values of the detuning parameter 0. From top 
to bottom = 1.3,1.3047,1.30478592,1.304788. I s = 0.9. 
Right: Sketch of the phase space for each parameter value. 
The thick line shows the trajectory of the LS in phase space. 



arbitrary precision from a semi-analytical stability anal- 
ysis of the unstable LS as in 8]. The lower-branch LS 
has one single positive eigenvalue Ai = 0.17713581. In 
Fig. [3| (right) we plot using crosses the period of the os- 
cillation LS as a function of ln(6 c — 6) obtained from 
numerical simulations of Eq. Q . Performing a linear fit- 
ting we obtain a slope 5.60, in excellent agreement with 
the theoretical prediction 1/Ai = 5.65, proving the exis- 
tence of a saddle- loop bifurcation for the oscillating LS. 
We should note that the theory was developed for planar 
bifurcations, therefore, as the phase space is a plane, the 
saddle has one unstable direction and one attracting di- 
rection pjj- The stable manifold of the unstable LS is, 
however, infinite dimensional. The success of the planar 
theory to describe our infinite dimensional system can be 
attributed to the fact that, somehow, the dynamics of the 
LS is effectively two-dimensional with a single unstable 
manifold and an effective stable manifold. 

In systems without spatial dependence it has been 
shown that an scenario composed by a saddle-loop bi- 
furcation and a stable fixed point leads to a excitability 
regime 0, Q, In our infinite-dimensional system 
LS does indeed show en excitable behavior: the homo- 
geneous solution is a globally attracting fixed point but 
localized disturbances (above the lower-branch LS) can 
send the system on a long excursion through phase space 
before returning to the fixed point. 

Fig. 01 shows the resulting trajectories of applying a 
perturbation in the direction of the unstable LS with 
three different intensities: one below the excitability 
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FIG. 3: a) Period of the limit cycle T as a function of the 
detuning for I s = 0.9. The vertical dashed line indicate the 
threshold of the saddle- loop bifurcation C = 1.30478592. b) 
Period T as a function of ln(# c — 6). Crosses correspond to 
numerical simulations while the solid line has a slope 1/Ai 
with Ai = 0.177 obtained from the linear stability analysis of 
the lower-branch LS. 
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threshold (dotted line), and two above: one very close to 
threshold (dashed line) and the other well above (solid 
line). In the first case the system relax exponentially to 
the homogeneous solution, while in the latter two cases 
it perform a long excursion before returning to the stable 
fixed point. In the case of a near threshold excitation the 
refractory period is appreciably longer due to the effect of 
the saddle. The spatial profile of the localized structure 
is shown in Fig. 0] The peak grows to a large value until 
the losses stop it. Then it decays exponentially until it 
disappears. A remnant wave is emitted out of the center 
dissipating the remaining energy. 

In absence of spatial degrees of freedom, Eq. JIJ does 
not present any kind of excitability. This behavior is 
strictly related to the dynamics of the 2D LS. Eq. JJ) is, 
in fact, a nonlinear Schrodinger equation (NLSE) with 
driving and damping. The phenomenon of self-focusing 
collapse |2£| of the 2D NLSE is behind the long excursion 
in phase space. When a localized perturbation concen- 
trates enough power, the self- focusing nonlinear mecha- 
nism induces a concentration of energy at that place. In 
the absence of losses a collapse, i.e. a divergence of the 
energy concentration, would take place in a finite time 
|2lj . However, the presence of losses prevents this col- 
lapse [2^ • The perturbation is finally dissipated and the 
system returns to the homogeneous solution. The same 
mechanism has also been used to explained the transport 
properties and dissipation rates in a wide class of turbu- 
lent flows |23|. The mechanism leading to excitability 
proposed in this letter is therefore not restricted to non- 
linear optics. It may have implications in plasma physics 
and hydrodynamics, where coherent structures may have 
similar features as LS. Furthermore, LS appearing on dif- 
ferent systems may undergo a Hopf bifurcation due to the 
nonlinear dynamics, even in the ID case, as for example 
in parametrically amplified optical systems j§| . 

In the limit of large detuning, the saddle-node, Hopf 
and saddle-loop bifurcation lines meet asymptotically, at 



FIG. 4: Evolution starting from the homogeneous solution 
(I s = 0.9) plus a localized perturbation of the form of the 
unstable LS multiplied by a factor a. Top panel shows the 
time evolution of the maximum intensity for a — 0.8 (dotted 
line); a = 1.01 (dashed line); a = 1.2 (solid line). The 3D 
plots show the transverse profile at different times for a — 
1.01. 



I s = as shown in Fig. ^ It is known that the intersec- 
tion of a saddle-node line with a Hopf line is a Takens- 
Bogdanov (TB) codimension-2 bifurcation point |24| iff 
there is a double zero eigenvalue (the imaginary part of 
the Hopf vanishes as approaching the intersection with 
the saddle- node). The unfolding around a TB point leads 
to a saddle- loop bifurcation line [24J. So, this unfolding 
fully explains the observed scenario, where once again 
our formally infinite-dimensional system appears to be 
perfectly described by a dynamical system in the plane. 

In Fig. we plot the two eigenvalues with largest real 
part of the LS for parameter values corresponding to 
three vertical cuts of Fig. QJl). Open symbols indicate 
eigenvalues with a non-zero imaginary part while filled 
symbols correspond to real eigenvalues. Where open 
symbols cross zero in Fig. (5Jl) corresponds to the Hopf 
bifurcation while filled symbols crossing zero indicates 
the saddle-node bifurcation. For the three plots we have 
taken as origin the saddle-node bifurcation. The impor- 
tant point is that in the line of the two complex conjugate 
eigenvalues responsible of the Hopf bifurcation there is a 
place where the imaginary part vanishes leading to two 
real eigenvalues, the largest of which is precisely the re- 
sponsible of the saddle-node bifurcation. As detuning in- 
creases the Hopf and saddle-node bifurcation points gets 
closer and closer but the structure of eigenvalues remains 
unchanged so that when the Hopf and saddle-node bifur- 
cation will finally meet the Hopf bifurcation will have 
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FIG. 5: Real (a) and imaginary (b) parts of the stable LS 
eigenvalues for three vertical cuts in Fig. [I] corresponding to 
detuning 0: 1.7 (□); 1.5 (A); 1.4 (0) versus the difference 
between I s and its value at the saddle- node bifurcation if(#). 



zero frequency, signaling a TB point. 

The TB point occurs asymptotically in the limit of 
large detuning and small pump Eq. This limit corre- 
sponds to the case in which Eq. becomes the conser- 
vative NLSE 25]. Details of the Hopf instability in this 
limit were studied in |l8|, where there is evidence of the 
double-zero bifurcation point, but the unfolding leading 
to the scenario presented here is not analyzed. 

In conclusion, we have analyzed an excitable regime 
associated to the existence of localized structures. We 
thus show that, in order to exhibit excitability, extended 
systems do not necessarily require to have local excitable 
behavior, instead such phenomenon can emerge due to 
the spatial dependence through the dynamics of a coher- 
ent (localized) structure. This open the possibility to ob- 
serve excitable behavior in a whole new class of systems 
where excitability was not thought to be present. Coher- 
ent resonance [2f| is also expected to be observed in the- 
ses systems upon application of localized disturbances of 
stochastic amplitude. Finally, localized structures have 
been shown to have a great potential as bits for optical 
memories 0, Q . This new excitable regimen opens also 
a new possibility for the use of localized structures as 
centers to process optical information in a similar way 
neurons do with electrical signals. 
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